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The quantum clock synchronization (QCS) algorithm proposed by Chuang (Phys. Rev. Lett, 85, 



<N , Abstract 

> 

00 

o 

vq 2006(2000)) has been implemented in a three qubit nuclear magnetic resonance quantum system. 
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The time difference between two separated clocks can be determined by measuring the output 



states. The experimental realization of the QCS algorithm also demonstrates an application of the 
quantum phase estimation. 
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I. INTRODUCTION 

The combination of quantum mechanics and computer science gives birth to quantum 
computer where quantum properties enable quantum computer to efficiently solve difficult 
problems in classical computer, for instance to factorize a large number using the Shor 
algorithm jl|. Marriage of quantum mechanics with other traditional science and technology 
has also produced fruitful results. Clock synchronization is one such example. The task of 
clock synchronization is to determine the time difference A between two spatially separated 
clocks. It is important both in practical application and in scientific research. Recently, 
some authors have applied quantum mechanical means to this classical problem and showed 
that quantum clock synchronization gains significant improvement compared to its classical 
counterpart j^, |3|, |4 1 • One of the quantum clock synchronization(QCS) algorithms is the one 
proposed by Chuang J2|, which uses the quantum phase estimation methodp, |6(. Chuang's 
QCS algorithm obtains n digits of accuracy in the time difference A while exchanging only 
the order of n number of qubits. This quantum algorithm gains an exponential improvement 
over classical algorithms which requires 0(2 n ) message exchanges. 

With the highest speed, photon is the natural choice for the practical implementation of 
the QCS algorithm. An alternative quantum system to implement the QCS algorithm is to 
use quantum spins in a magnetic field as suggested in Ref . J2| • Nuclear magnetic resonance 
(NMR) has been widely used in various fields, and it has also become an important arena 
to demonstrate quantum algorithms. Many algorithms, such as the Grover algorithm, the 
quantum Fourier transform, and Shor's quantum factoring algorithm, have been demon- 



strated in NMR quantum systems [7|, l8ll9lllC. 
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ll( . Some quantum communication protocols, 



such as the quantum teleportation, quantum dense coding have also been demonstrated 
in NMR quantum systems[12|, |l3|, Il4j . It is thus a good system to demonstrate the QCS 
algorithm. 

In this paper, we report the result of an implementation of the QCS algorithm in a 
three qubit NMR quantum computer. The basis of the algorithm is the quantum Fourier 
transform (QFT) which has been applied in experiments previously JQL I la. Il6l|. The QCS 
algorithm requires an pure state as the initial state. To use an NMR ensemble for quantum 
computation, preparation of the effective-pure state is one necessary step |l7|, ll8| . Temporal 
averaging and spatial averaging are two main practical methods to prepare effective-pure 
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2l|. In our experiments, we chose the spatial averaging method to prepare 
the effective-pure state. Using this method, the algorithm is implemented through only 
one experiment. Compared with the temporal averaging, the spatial ave rag i ng can shorten 
experiment time, and it has been applied in various experiments 22J, |23|, |2J, |25|, |26|, I27ll28j|. 
This paper is organized as follows. After this brief introduction, we give the outline of the 
QCS algorithm, and in particular we give the explict quantum circuit of the QCS algorithm 
in a three qubit system in section |H] In section IHIl we give the details of the experimental 
demonstration of the QCS algorithm, including the pulse sequence and the results of the 
experiment. In section HV| we give a brief summary. 

II. THE QCS ALGORITHM AND THE SIMPLIFIED CIRCUIT FOR A 3-QUBIT 
NMR QUANTUM SYSTEM 

A brief summary of the QCS algorithm is given below, and the details can be found 
in Ref. |2j. In the QCS algorithm, Alice possesses m + 1 qubits. The first m qubits are 
the working qubits and are retained at Alice's site, and the extra qubit is an ancilla qubit 
for implementing the T-unitary operation. The procedure of the QCS is as follows. Alice 
performs the Hadamard operation on each of the first m qubits. This prepares the state of 
the m + 1 qubit system to 

i<M = i 2 i:V>io>. 

V z fc=0 

Then the T operation, which inscribe the time difference between the two clocks into the 
many-body quantum state, is implemented. At this stage the state of m + 1 qubit system is 

|0 2 ) = -L 2 ^ 1 e 2 ^ A |A ; )|O). 
V ^ fe=o 

Then Alice applies an inverse quantum Fourier transform F~ l on the first m qubits, trans- 
forming their state to 

2 m -12 m -l 2 m — 1 

2 m 

The state is peaked at j = 2 m ujA. If 2 m u;A is an integer, the equality is exact. Measuring 
the qubit values of the first m qubits gives us the value of j which in turn determines the 
time difference A. In practice, the equality can not hold exactly, and it gives uA about m 
bits of accuracy with a high probability. 
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The central ingredients of the QCS algorithm are the T operation and the quantum 
Fourier transform. The T operation can be implemented by performing the following oper- 
ation on each of the m qubits: 1) Alice makes CNOT operation on the Z-th and the ancilla 
qubit which is the target state ; 2) The ticking qubit handshake protocol TQH(tc2 1 ~ 1 uj, \ip)) 
is performed, transforming the state of the Z-th and ancilla qubit system into 

V * 
3) Alice performs another CNOT gate to the Z-th and the ancila qubit, so that the state of 
the Z-th and the ancilla qubit becomes 
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2WA/2 j2 i^)e 2i * A |0>. (3) 
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4) After all the m qubit have gone through the previous operations, the overall operation is 
to transform the state \k')\Q) = \k ki ■ ■ -fc m _i)|0) into \k')e 27TiujA{ ^ 2 ' kl) \0), because 

T\k')\0) = [T \k )T 1 \k 1 )---T m _ 1 \k m _ 1 )]\0) = \k')e 2 ™ A ^i 2l ^\0). (4) 

Using the SWAP operation, the state \k') is transformed into \k) = \k m -i ■ ■ ■ kik ), and Eq. 
(J3J) becomes 

\k)e 2mkulA \0}, (5) 

where J2i^ l h = k. 

The quantum network shown in Fig. ^ implements the QCS algorithm in a three qubit 
system. The three lines denote the three qubits respectively. |0) denotes the spin up state. 
H denotes the Hadamard transform. The effect of the TQH is to introduce a phase to 
the two different quantum states of the ancilla qubit, namely for the |0) state with phase 
e -2WA/2 and for gtate [i) w ith phase e 2 ' niujA/2 } hence the state ^(|00) + |11)) of the Z-th 
qubit and the ancilla qubit system changes to 

1 ^-2'^A/2| 00) + e 2WA/2| n ^ (6) 

In NMR, this is equivalent to a rotation about the z-axis. For the first qubit, the rotation 
is R z {(pk) — e lipkIz , and for the second qubit the rotation is R z (2(pk) = R z (ipk)R z (Vk)- We 
have written — uAir as if. Because 2 2 t<;A has to be integer, ranging from 1 to 2 m — 1 = 3, so 
ip can take the following values (pk = —kn/2, where k — 0, 1, 2, 3. When m is large, k takes 



the value from to 2 m — 1 and the measurement of this k value gives the value of a; A, and 
hence the time difference A. 

For the quantum Fourier transform, some simplification is possible. We use IZ*J to 
denote the controlled phase shift operation applied to the subsystem constructed by qubit 
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1 and 2. L^J is explicitly 
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with the basis orders in |00), |01), |10), |11). The network outlined by the dashed rectangle in 
FigQ implements the inverse of the quantum fourier transform, where the SWAP operation 
has been counteracted by another one in the network. The inverse of the QFT F' 1 can be 
written as F /_1 SWAP|9|. The effect of this network is to make the following transformation 
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The bit values of the first two qubits are the desired output. The network shown in Fig. Q 
transforms |000) to |000), |010), 1 100) , and 1 110) , corresponding to ifk = 0, —n/2, — ir, and 
— 37r/2, respectively, and uA to be 0, 1/4, 2/4 and 3/4 respectively. By measuring qubit 1 
and 2, one obtains the concrete value of p and hence determines the time difference between 
the two clocks. 

III. IMPLEMENTATION IN A 3-QUBIT NMR QUANTUM SYSTEM 

The experiment uses a sample of Carbon- 13 labelled trichloroethylene (TCE) dissolved 
in d-chloroform. Data are taken at room temperature with a Bruker DRX 500 MHz spec- 
trometer. 1 H is denoted as qubit 3, the 13 C directly connecting to 1 H is denoted as qubit 2, 
and the other 13 C is denoted as qubit 1. The three qubits are denoted as CI, C2 and H3. 
By setting h = 1, the Hamitonian of the three-qubit system is |29] 
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P z (j = 1, 2, 3) are the matrices for ^-component of the angular momentum of the spins. V\, 
z/ 2 , V3 are the resonance frequencies of CI, C2 and H3, and v\ = v-i + 904.4Hz. The coupling 
constants are measured to be J i2 = 103.1 Hz, J 23 = 203.8 Hz, and J i3 = 9.16 Hz. The 
coupled-spin evolution between two spins is denoted as 

\r\. x = e -^J 3 iriU' % ( 10 ) 

where I = 1,2, 3, and j 7^ I. [r]ji can be realized by averaging the coupling constants other 
than Jji to zero 30]. For example, [r]i 3 is realized by the pulse sequence shown in Fig. El 
(a). The chemical shift evolution of C2 is realized by the pulse sequence shown in Fig. 
121 (b). R ^(tt) can be realized by choosing the proper evolution time, and the transmitter 
frequency J3l|. The it pulses for C2 are chosen as RE-BURP pulses to excite the multiplet 
of C2 uniformly 23- 

n 

The initial effective-pure state |000) is prepared by spatial averaging 20] . The following 
radio-frequency (rf ) pulse and gradient pulse sequence 



[tt/4] 1 / - [1/2J 12 ] 12 - [-57T/6];- 2 - [a]l - [grad], - [tt/4]J - [9/2J, 



23J23 — 

[l/2Ji 3 ]i3 - [tt/4]J - [grad], - [tt/4]J - [9/4J 23 ] 23 - [l/4J 13 ]i 3 - [ix/Af x - [grad],,(ll) 
transforms the system from the equilibrium state 

p eq = lc {ll + ID + JhI 3 z , (12) 

to 

p = ll/2 + ll/2 + ll/2 + 1\I 2 Z + l 2 z ll + l\l\ + 2l\l\ll (13) 

where an overall phase factor has been ignored J20j, |2l|, |22|, |3J|. [7r/4]^.' 2 denotes the 7r/4 
pulse exciting CI and C2 simultaneously along x-axis. [vr/4]^ denotes the spin-selective 
7r/4 pulse for X H along y-axis. jc an d Jh denotes the gyromagnetic ratio of 13 C and X H. 
a = arccos(— 7cV / 6/7h)- Po is equivalent to |000). We find that the compound operation 

CNOT 13 i^O fc )CNOT 13 = [^]ia, (14) 
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CNOT 23 i^(2^)CNOT 23 = [-^]23, (15) 



realize the network much easier. The Hadamard transform simultaneously applied to CI 
and C2, is realized by pulse sequence [— vr/2]^' 2 — [ir]],' 2 where the number in the superscript 
refers to the qubit number. The Hadamard transform for C2 in the inverse of QFT, denoted 
by H 2 , is realized by [tt/4]^ 2 ' 3 - R 2 z (tt) - [-tt/4]^ 2 ' 3 , and H 1 = H^H 2 , noting H' 1 = H. 
I"*/ 2 is realized by [1/4J 12 ] - [-vr/2]^ - [vr/4]^ 2 - [tf/^QQ. 

The experimental results are represented as the density matrices obtained by the state 
tomography technique |l9j, |35j, [36] , where the spin-selective readout pulses for C2, denoted 
as [vr/2] 2 and [7r/2] 2 are realized by pulse sequences [vr/2]^' 2 — R 2 z (n/2) — \—7t/ 2}}> 2 and 



-7r/2]J.' 2 — R 2 z {ti/2) — [7r/2]^' 2 , respectively, in order to freeze the motion of CI 



jli 2 and 



R 2 (tt/2) is realized by the pulse sequence shown in Fig. |21^b). Fig. |S] shows the experimen- 
tally measured density matrix when the system lies in effective- pure state |000) prepared by 
the pulse sequence (jllj) . In the generated density matrix, the desired element, which is the 
only nonzero element in theory, is measured to be 14.2 (in arbitrary units). The amplitudes 
of the other elements, which are zero in theory, are less than 1.8. 

The QCS algorithm starts with effective- pure state|000). When <p — 0, (f\ — —n/2, 
(f2 = —n, and if 3 = — 37r/2, the network shown in Fig. Q transforms |000) to |000), |010), 
1 100) , and 1 110) , respectively, corresponding to the four different time differences A = 0, 
A = l/4a;, A = l/2u>, and A = 3 /4a;. Figs. 0] (a-d) show the experimentally measured 
density matrices of the three- qubit system after the completion of the QCS algorithm, 
corresponding to (po-^s, respectively. The fidelity of the transformation is described by [2| 
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Pinitiai is the initial density matrix shown in Fig. El Ptheory = Up ini ti a iU^ , where U denotes the 
theoretical transformation to implement the QCS algorithm. p exp denotes the experimentally 
measured density matrix shown in Figs. EJ The fidelities corresponding to Figs. 0] (a- 
d) are 90.7%, 77.4%, 75.2%, and 77.0%, respectively. The errors mainly result from the 
imperfection of the pulses, the inhomogeneity in the magnetic field and the decoherence 
time limit. 



IV. CONCLUSION 

We have implemented the quantum clock synchronization algorithm in a three- qubit 
NMR quantum computer. Using spatial averaging, the algorithm is implemented through 
one experimental run which is much shorter than the temporal averaging method. For small 
qubit system such as three qubit used in this experiment, the level of accuracy as good as 
the temporal averaging method. Compared with temporal averaging, process of experiments 
is simplified greatly. The time difference can be can be read out through the output state. 
In the experiments, we have exploited the long range coupling between non-adjacent nuclear 
spins, and it is found that it works well. In NMR samples, long range coupling is a precious 
resource, and use of this long range coupling in experimental realization should be made 
as much as possible. On the other hand, through optimizing network, the experimental 
operation difficulty can be reduced. In this experiment, using a simplified network where 
some redundant operations have been got rid of, we have shortened the total time consumed 
by the whole algorithm, and hence also reduces the effect of decoherence consequently. 
Although we have demonstrated the algorithm with only three qubit, the techniques can be 
generalized to more qubit system. It should be pointed out QCS is one promising area of 
quantum information technology that may be implemented in the future because it requires 
less number of qubit than other quantum algorithms. For instance to achieve an accuracy of 
lOOps, the number of qubits required is only 34. Of course, before QCS can be realistically 
used in practice, there should be big advancement in the technology, for instance in the 
accuracy of quantum gate operations, and the length of decoherence time. 
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FIG. 1: The quantum network to implement the QCS algorithm. The three lines denote the three 

— 7r/2 

qubits. H denotes the Hadmamard transform, I, n / denotes the controlled phase shift operation, 
and R z ((fk) = e i<Pk , where ifk = —kir/2 (k = 0,1,2,3). The network outlined by the dashed 
rectangle implements the inverse of quantum Fourier transform without the SWAP operation. 
Time goes from left to right. |x)i|y)2 is the output state, which can be obtained through measuring 
qubits 1 and 2. 
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FIG. 2: Pulse sequence to realize [r]i3 (a) and the chemical shift evolution R^,(2Ti58t) (b), where 
S is the offset between the chemical shift of C2 and the transmitter frequency. Rectangles denote 
7r pulses. The narrow ones are so short that the widths can be ignored. The wide ones, however, 
are long pulses of which widths can not be ignored. 
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FIG. 3: The density matrix (in arbitrary units) reconstructed using the state tomography technique, 
when the system is prepared in effective- pure state |000). 
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FIG. 4: The density matrices of the three- qubit system, after the implementation of the QCS 
algorithm. Only the real components are plotted. The imaginary portions, which are theoretically 
zero, are found to contribute less than 10% to the experimental results. Figs, (a-d) correspond 
to the four effective- pure states |000), |010), 1 100) , and 1 110) , and the four effective- pure states 
correspond to the four different time differences A = 0, A = 1/4cj, A = l/2u, and A = 3/4u, 
respectively. 
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